Spatial structure, chemotaxis and quorum sensing shape bacterial biomass accumulation in complex porous media.

Biological tissues, sediments, or engineered systems are spatially structured media with a tortuous and porous structure that host the flow of fluids. Such complex environments can influence the spatial and temporal colonization patterns of bacteria by controlling the transport of individual bacterial cells, the availability of resources, and the distribution of chemical signals for communication. Yet, due to the multi-scale structure of these complex systems, it is hard to assess how different biotic and abiotic properties work together to control the accumulation of bacterial biomass. Here, we explore how flow-mediated interactions allow the gut commensal Escherichia coli to colonize a porous structure that is composed of heterogenous dead-end pores (DEPs) and connecting percolating channels, i.e. transmitting pores (TPs), mimicking the structured surface of mammalian guts. We find that in presence of flow, gradients of the quorum sensing (QS) signaling molecule autoinducer-2 (AI-2) promote E. coli chemotactic accumulation in the DEPs. In this crowded environment, the combination of growth and cell-to-cell collision favors the development of suspended bacterial aggregates. This results in hot-spots of resource consumption, which, upon resource limitation, triggers the mechanical evasion of biomass from nutrients and oxygen depleted DEPs. Our findings demonstrate that microscale medium structure and complex flow coupled with bacterial quorum sensing and chemotaxis control the heterogenous accumulation of bacterial biomass in a spatially structured environment, such as villi and crypts in the gut or in tortuous pores within soil and filters.

A Method to Measure the Diffusion Coefficient in Liquids.

Molecular diffusion in liquids is a key process in numerous systems: it is often the reaction rate limiting factor in biological or chemical reaction. Molecular diffusion has been recognized as the ultimate mechanism by which substances concentration get homogenized and, thus, their mixing and dilution occur. Here, we propose a novel method to directly measure the diffusion coefficient D of solutes or suspensions in liquids. Differently from current methods, as Dynamic Light Scattering or Fluorescent Correlation Spectroscopy, our method does not rely on previous knowledge on the fluid or tracer properties, but it is based on directly measuring the concentration spatial profile of a considered tracer with optical techniques within a diffusion chamber. We test this novel method on a sample of mono-dispersed suspension of spherical colloids for which an estimate for D can be made based on Einstein-Stokes relation. We, then, use this technique to measure the diffusion coefficient of a non-spherical tracer. We further quantify mixing of the considered tracers in the confined domain of the diffusion chamber: we show that, since diffusion-limited mixing (quantified in terms of the dilution index) in a confined space happens faster than un-confined domain, the finite size of the diffusion chamber must be taken into account to properly estimate D and the tracer mixing degree.

Electrical Signatures of Diffusion-Limited Mixing: Insights from a Milli-fluidic Tracer Experiment.

We investigate how diffusion-limited mixing of a layered solute concentration distribution within a porous medium impacts bulk electrical conductivity. To do so, we perform a milli-fluidic tracer test by injecting a fluorescent and electrically conductive tracer in a quasi two-dimensional (2D) water-saturated porous medium. High resolution optical- and geoelectrical monitoring of the tracer is achieved by using a fluorimetry technique and equipping the flow cell with a resistivity meter, respectively. We find that optical and geoelectrical outputs can be related by a temporal re-scaling that accounts for the different diffusion rates of the optical and electrical tracers. Mixing-driven perturbations of the electrical equipotential field lines cause apparent electrical conductivity time-series, measured perpendicularly to the layering, to peak at times that are in agreement with the diffusion transport time-scale associated with the layer width. Numerical simulations highlight high sensitivity of such electrical data to the layers' degree of mixing and their distance to the injection electrodes. Furthermore, the electrical data correlate well with time-series of two commonly used solute mixing descriptors: the concentration variance and the scalar dissipation rate.

Morphology and Size of Bacterial Colonies Control Anoxic Microenvironment Formation in Porous Media.

Bacterial metabolisms using electron acceptors other than oxygen (e.g., methanogenesis and fermentation) largely contribute to element cycling and natural contaminant attenuation/mobilization, even in well-oxygenated porous environments, such as shallow aquifers. This paradox is commonly explained by the occurrence of small-scale anoxic microenvironments generated by the coupling of bacterial respiration and dissolved oxygen (O2) transport by pore water. Such microenvironments allow facultative anaerobic bacteria to proliferate in oxic environments. Microenvironment dynamics are still poorly understood due to the challenge of directly observing biomass and O2 distributions at the microscale within an opaque sediment matrix. To overcome these limitations, we integrated a microfluidic device with transparent O2 planar optical sensors to measure the temporal behavior of dissolved O2 concentrations and biomass distributions with time-lapse videomicroscopy. Our results reveal that bacterial colony morphology, which is highly variable in flowing porous systems, controls the formation of anoxic microenvironments. We rationalize our observations through a colony-scale Damköhler number comparing dissolved O2 diffusion and a bacterial O2 uptake rate. Our Damköhler number enables us to predict the pore space fraction occupied by anoxic microenvironments in our system for a given bacterial organization.

Microfluidic study in a meter-long reactive path reveals how the medium's structural heterogeneity shapes MICP-induced biocementation.

Microbially induced calcium carbonate (CaCO3) precipitation (MICP) is one of the major sustainable alternatives to the artificial cementation of granular media. MICP consists of injecting the soil with bacterial- and calcium-rich solutions sequentially to form calcite bonds among the soil particles that improve the strength and stiffness of soils. The performance of MICP is governed by the underlying microscale processes of bacterial growth, reactive transport of solutes, reaction rates, crystal nucleation and growth. However, the impact of pore-scale heterogeneity on these processes during MICP is not well understood. This paper sheds light on the effect of pore-scale heterogeneity on the spatiotemporal evolution of MICP, overall chemical reaction efficiency and permeability evolution by combining two meter-long microfluidic devices of identical dimensions and porosity with homogeneous and heterogeneous porous networks and real-time monitoring. The two chips received, in triplicate, MICP treatment with an imposed flow and the same initial conditions, while the inlet and outlet pressures were periodically monitored. This paper proposes a comprehensive workflow destined to detect bacteria and crystals from time-lapse microscopy data at multiple positions along a microfluidic replica of porous media treated with MICP. CaCO3 crystals were formed 1 h after the introduction of the cementation solution (CS), and crystal growth was completed 12 h later. The average crystal growth rate was overall higher in the heterogeneous porous medium, while it became slower after the first 3 h of cementation injection. It was found that the average chemical reaction efficiency presented a peak of 34% at the middle of the chip and remained above 20% before the last 90 mm of the reactive path for the heterogeneous porous network. The homogeneous porous medium presented an overall lower average reaction efficiency, which peaked at 27% 420 mm downstream of the inlet and remained lower than 12% for the rest of the microfluidic channel. These different trends of chemical efficiency in the two networks are due to a higher number of crystals of higher average diameter in the heterogeneous medium than in the homogeneous porous medium. In the interval between 480 and 900 mm, the number of crystals in the heterogeneous porous medium is more than double the number of crystals in the homogeneous porous medium. The average diameters of the crystals were 23-46 µm in the heterogeneous porous medium, compared to 17-40 µm in the homogeneous porous medium across the whole chip. The permeability of the heterogeneous porous medium was more affected than that of the homogeneous system, while the pressure sensors effectively captured a higher decrease in the permeability during the first two hours when crystals were formed and a less prominent decrease during the subsequent seeded growth of the existing crystals, as well as the nucleation and growth of new crystals.

Structure induced laminar vortices control anomalous dispersion in porous media.

Natural porous systems, such as soil, membranes, and biological tissues comprise disordered structures characterized by dead-end pores connected to a network of percolating channels. The release and dispersion of particles, solutes, and microorganisms from such features is key for a broad range of environmental and medical applications including soil remediation, filtration and drug delivery. Yet, owing to the stagnant and opaque nature of these disordered systems, the role of microscopic structure and flow on the dispersion of particles and solutes remains poorly understood. Here, we use a microfluidic model system that features a pore structure characterized by distributed dead-ends to determine how particles are transported, retained and dispersed. We observe strong tailing of arrival time distributions at the outlet of the medium characterized by power-law decay with an exponent of 2/3. Using numerical simulations and an analytical model, we link this behavior to particles initially located within dead-end pores, and explain the tailing exponent with a hopping across and rolling along the streamlines of vortices within dead-end pores. We quantify such anomalous dispersal by a stochastic model that predicts the full evolution of arrival times. Our results demonstrate how microscopic flow structures can impact macroscopic particle transport.

Morphogenesis and oxygen dynamics in phototrophic biofilms growing across a gradient of hydraulic conditions.

Biofilms are surface-attached and matrix-enclosed microbial communities that dominate microbial life in numerous ecosystems. Using flumes and automated optical coherence tomography, we studied the morphogenesis of phototrophic biofilms along a gradient of hydraulic conditions. Compact and coalescent biofilms formed under elevated bed shear stress, whereas protruding clusters separated by troughs formed under reduced shear stress. This morphological differentiation did not linearly follow the hydraulic gradient, but a break point emerged around a shear stress of ~0.08 Pa. While community composition did not differ between high and low shear environments, our results suggest that the morphological differentiation was linked to biomass displacement and reciprocal interactions between the biofilm structure and hydraulics. Mapping oxygen concentrations within and around biofilm structures, we provide empirical evidence for biofilm-induced alterations of oxygen mass transfer. Our findings suggest that architectural plasticity, efficient mass transfer, and resistance to shear stress contribute to the success of phototrophic biofilms.

Combining Fluidic Devices with Microscopy and Flow Cytometry to Study Microbial Transport in Porous Media Across Spatial Scales.

Understanding the transport, dispersion and deposition of microorganisms in porous media is a complex scientific task comprising topics as diverse as hydrodynamics, ecology and environmental engineering. Modeling bacterial transport in porous environments at different spatial scales is critical to better predict the consequences of bacterial transport, yet current models often fail to up-scale from laboratory to field conditions. Here, we introduce experimental tools to study bacterial transport in porous media at two spatial scales. The aim of these tools is to obtain macroscopic observables (such as breakthrough curves or deposition profiles) of bacteria injected into transparent porous matrices. At the small scale (10-1000 µm), microfluidic devices are combined with optical video-microscopy and image processing to obtain breakthrough curves and, at the same time, to track individual bacterial cells at the pore scale. At larger scale, flow cytometry is combined with a self-made robotic dispenser to obtain breakthrough curves. We illustrate the utility of these tools to better understand how bacteria are transported in complex porous media such as the hyporheic zone of streams. As these tools provide simultaneous measurements across scales, they pave the way for mechanism-based models, critically important for upscaling. Application of these tools may not only contribute to the development of novel bioremediation applications but also shed new light on the ecological strategies of microorganisms colonizing porous substrates.

Hyperuniform Monocrystalline Structures by Spinodal Solid-State Dewetting.

Materials featuring anomalous suppression of density fluctuations over large length scales are emerging systems known as disordered hyperuniform. The underlying hidden order renders them appealing for several applications, such as light management and topologically protected electronic states. These applications require scalable fabrication, which is hard to achieve with available top-down approaches. Theoretically, it is known that spinodal decomposition can lead to disordered hyperuniform architectures. Spontaneous formation of stable patterns could thus be a viable path for the bottom-up fabrication of these materials. Here, we show that monocrystalline semiconductor-based structures, in particular Si_{1-x}Ge_{x} layers deposited on silicon-on-insulator substrates, can undergo spinodal solid-state dewetting featuring correlated disorder with an effective hyperuniform character. Nano- to micrometric sized structures targeting specific morphologies and hyperuniform character can be obtained, proving the generality of the approach and paving the way for technological applications of disordered hyperuniform metamaterials. Phase-field simulations explain the underlying nonlinear dynamics and the physical origin of the emerging patterns.

Trait-specific dispersal of bacteria in heterogeneous porous environments: from pore to porous medium scale.

The dispersal of organisms controls the structure and dynamics of populations and communities, and can regulate ecosystem functioning. Predicting dispersal patterns across scales is important to understand microbial life in heterogeneous porous environments such as soils and sediments. We developed a multi-scale approach, combining experiments with microfluidic devices and time-lapse microscopy to track individual bacterial trajectories and measure the overall breakthrough curves and bacterial deposition profiles: we, then, linked the two scales with a novel stochastic model. We show that motile cells of Pseudomonas putida disperse more efficiently than non-motile mutants through a designed heterogeneous porous system. Motile cells can evade flow-imposed trajectories, enabling them to explore larger pore areas than non-motile cells. While transported cells exhibited a rotation in response to hydrodynamic shear, motile cells were less susceptible to the torque, maintaining their body oriented towards the flow direction and thus changing the population velocity distribution with a significant impact on the overall transport properties. We also found, in a separate set of experiments, that if the suspension flows through a porous system already colonized by a biofilm, P. putida cells are channelled into preferential flow paths and the cell attachment rate is increased. These two effects were more pronounced for non-motile than for motile cells. Our findings suggest that motility coupled with heterogeneous flows can be beneficial to motile bacteria in confined environments as it enables them to actively explore the space for resources or evade regions with unfavourable conditions. Our study also underlines the benefit of a multi-scale approach to the study of bacterial dispersal in porous systems.

Unraveling the biophysical underpinnings to the success of multispecies biofilms in porous environments.

Biofilms regulate critical processes in porous ecosystems. However, the biophysical underpinnings of the ecological success of these biofilms are poorly understood. Combining experiments with fluidic devices, sequencing and modeling, we reveal that architectural plasticity enhances space exploitation by multispecies biofilms in porous environments. Biofilms consistently differentiated into an annular base biofilm coating the grains and into streamers protruding from the grains into the pore space. Although different flow-related processes governed the differentiation of these architectures, both BB and streamers were composed of similar bacterial assemblages. This is evidence for architectural plasticity. Architectural plasticity allowed for complementary use of the space provided by the grain-pore complexes, which increased biofilm carrying capacity at the larger scale of the porous system. This increase comes potentially at the cost of a tradeoff. Contrasting time scales of oxygen replenishment and consumption, we show that streamers locally inhibit the growth of the BB downstream from the grains. Our study provides first insights into the biophysical underpinnings to the success of multispecies biofilms in porous environments.

Interface evolution during radial miscible viscous fingering.

We study experimentally the miscible radial displacement of a more viscous fluid by a less viscous one in a horizontal Hele-Shaw cell. For the range of tested injection rates and viscosity ratios we observe two regimes for the evolution of the fluid-fluid interface. At early times the interface length increases linearly with time, which is typical of the Saffman-Taylor instability for this radial configuration. However, as time increases, the interface growth slows down and scales as ∼t(1/2), as one expects in a stable displacement, indicating that the overall flow instability has shut down. Surprisingly, the crossover time between these two regimes decreases with increasing injection rate. We propose a theoretical model that is consistent with our experimental results, explains the origin of this second regime, and predicts the scaling of the crossover time with injection rate and the mobility ratio. The key determinant of the observed scalings is the competition between advection and diffusion time scales at the displacement front, suggesting that our analysis can be applied to other interfacial-evolution problems such as the Rayleigh-Bénard-Darcy instability.

Flow intermittency, dispersion, and correlated continuous time random walks in porous media.

We study the intermittency of fluid velocities in porous media and its relation to anomalous dispersion. Lagrangian velocities measured at equidistant points along streamlines are shown to form a spatial Markov process. As a consequence of this remarkable property, the dispersion of fluid particles can be described by a continuous time random walk with correlated temporal increments. This new dynamical picture of intermittency provides a direct link between the microscale flow, its intermittent properties, and non-Fickian dispersion.

Anomalous kinetics in diffusion limited reactions linked to non-Gaussian concentration probability distribution function.

We investigate anomalous reaction kinetics related to segregation in the one-dimensional reaction-diffusion system A + B → C. It is well known that spatial fluctuations in the species concentrations cause a breakdown of the mean-field behavior at low concentration values. The scaling of the average concentration with time changes from the mean-field t(-1) to the anomalous t(-1/4) behavior. Using a stochastic modeling approach, the reaction-diffusion system can be fully characterized by the multi-point probability distribution function (PDF) of the species concentrations. Its evolution is governed by a Fokker-Planck equation with moving boundaries, which are determined by the positivity of the species concentrations. The concentration PDF is in general non-Gaussian. As long as the concentration fluctuations are small compared to the mean, the PDF can be approximated by a Gaussian distribution. This behavior breaks down in the fluctuation dominated regime, for which anomalous reaction kinetics are observed. We show that the transition from mean field to anomalous reaction kinetics is intimately linked to the evolution of the concentration PDF from a Gaussian to non-Gaussian shape. This establishes a direct relationship between anomalous reaction kinetics, incomplete mixing and the non-Gaussian nature of the concentration PDF.

Spatial model of autocatalytic reactions.

Biological cells with all of their surface structure and complex interior stripped away are essentially vesicles--membranes composed of lipid bilayers which form closed sacs. Vesicles are thought to be relevant as models of primitive protocells, and they could have provided the ideal environment for prebiotic reactions to occur. In this paper, we investigate the stochastic dynamics of a set of autocatalytic reactions, within a spatially bounded domain, so as to mimic a primordial cell. The discreteness of the constituents of the autocatalytic reactions gives rise to large sustained oscillations even when the number of constituents is quite large. These oscillations are spatiotemporal in nature, unlike those found in previous studies, which consisted only of temporal oscillations. We speculate that these oscillations may have a role in seeding membrane instabilities which lead to vesicle division. In this way synchronization could be achieved between protocell growth and the reproduction rate of the constituents (the protogenetic material) in simple protocells.